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Abstract 

We use the properties of Hermite and Kampe de Feriet polynomials to get closed forms for the 
repeated derivatives of functions whose argument is a quadratic or higher-order polynomial. The 
results we obtain are extended to product of functions of the above argument, thus giving rise 
to expressions which can formally be interpreted as generalizations of the familiar Leibniz rule. 
Finally, examples of practical interest are discussed. 
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I. INTRODUCTION 

Formula (1.1.1.1) of Ref. [ij] refers to repeated derivatives of functions whose argument 
is a quadratic polynomial, i.e. {D^ = d/d^) 

[n/2] \n-2k 

4"[/(-^)l="!g^f^;^«:r'=/(.^). (1.1) 

Albeit elementary, and a particular case of the Faa di Bruno formula some aspects of 
Eq. (II. ip deserve to be studied as they may lead to novel results. 

Eq. f ll.ip implicitly assumes that f{x) is at least a n-times different iable function. More- 
over, we make the hypothesis that it can be written as 

/(^)=E^/- = «'^/0, rf0 = fn, (1.2) 

where /o denotes a kind of "vacuum" state on which the repeated action of the umbral 
operator (p generates a discrete sequence of functions. In what follows we will consider only 
functions that admit the expansion ( 11.21) . i.e., that are analytical over the entire complex 
plane. 

As a. example, we note *a. .he T,.om, fu„c«o„ of o.der .ero B can be wnUen as 

^o(^) = E7^ = ^"^'/o, (1.3) 

with 

^'f' = WTTy A=p^ = i. (1.4) 

The operator may also be raised to any real (not necessarily integer) power, so that the 
Tricomi function of order a can be written as 

The obvious advantage of the previous umbral re-shaping of the function f{x) is the possi- 
bility of exploiting the wealth of the properties of the exponential function, which will be 
assumed to be still valid, since we will treat the operator as an ordinary constant [sl. 
As a consequence of Eq. ( 11.21) , we can write (we omit the vacuum /o for brevity) 

^,%)/[^?(^)] = 0"e^(-)^. (1.6) 



By taking into account the identity 4| 

D^e'^"' = i7^')(2ax,a)e'^"', (1.7) 

where 

[n/2] 

{n-2k)\k\ 



A;=0 



is the two-variable Hermite-Kampe de Feriet polynomials 



blf{x') = ble-'^ = iff (2x0,0)e^ 



Q3, 



one obtains 



^1 y e-'^ (1.9) 

f^^ {n - 2k)\ k\ ^ I ' ^ ^ 

from which, by using Eq. (II. 6^ . we recover Eq. 01. ip (see also Appendix A). 

II. AN EXTENSION OF LEIBNIZ FORMULA 

In this section we will draw further consequences from the previous formalism. In par- 
ticular, we will consider the case in which the function f{x) is the product of two functions, 
obtaining a closed formula which will be recognized as a generalized version of the Leibniz 
rule. 

We start by considering the following function 

f{x')=g{x')h{x') = (e^'^'go) (e^'^/^o) , (H-l) 
i.e., since the operators 7 and 17 commute (we omit the vacuum terms) 

/(x^) =e"'(^+^^ (II.2) 
and, according to the Eq. ( 11. 9^ . we find 

b^lfix'')] = £);^e"'(^+^^) = iff (2(7 + r})a;,7 + r}) e"'(^+^) 

"'E^i^R + (.1.3) 



(n-2k)lkl ^\ 3 
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By noting that 



we obtain 



[n/2] 



,n-2fc 



n—k 



n-2k)\k\ 

k=0 ^ ' j=0 



n — k 
J 



4™ 9ix'] 



(11.4) 



(11.5) 



which is manifestly a generahzation of the ordinary Leibniz rule for the nth derivative of 
the product of two functions. The use of the addition formula [5| 

Hi'\x^+X2,yi + y2) = Q H;^lkix^,y^) Hi'\x2,y2) (II.6) 

fc=0 ^ ^ 

in Eq. fllL3|) yields 

n 

D:[g{x')h{x')] = j2 



k=0 



ifi5,(2x7,7)e^'^ Hlf>{2xfj,fj)e^^^ 



(2), 



(11.7) 



which makes the analogy more transparent. From this example, we conclude that, for the 
type of function we are dealing with, the operator 

(2)f(«) = if(2)(2 7x,7) (II.8) 

is a kind of multiple derivative operator. For more detail on this point see Appendix B. 

A fairly interesting application of Eq. (III.3P is in the calculation of nth derivative of the 
product of cylindrical Bessel functions, that, in our present formalism, writes 



JJx) 



X 



e-(-/^)^^jo, 



^ 3o= 3r. 



r(n + l)' 



One obtains 



[n/2] 



D:[Ji{x)] = {-irn\ J2 



'-2xy 



n—k 



n — k 



m 



(11.9) 



(11.10) 



{n-2k)\k\ 

Obviously, the formalism can easily applied also to other products of functions. 

We close this section, by considering n-th derivative of the function /(ax^ + hx). From 
Eq. (11.2p . one has 

bl[f{ax'^ + hx)] =i)!^e("^'+^^-)*/o (11.11) 

and the use of the identity 

jjn^ax^+hx ^ ij^2)^2ax + 6,a)e"^'+^^ (11.12) 

leads to 

(9.nT A- h\r^-2k ak ^ , 

(11.13) 



D:[fiax^ + bx)] =nlY: ^ ^ L ^^^.(-^ + 



fc=0 



{n-2k)\k\ 
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III. HIGHER-ORDER HERMITE POLYNOMIALS AND REPEATED DERIVA- 
TIVES 



We have already remarked that the examples we are considering here are the particular 
cases of the Faa di Bruno formula, which had different formulations in the course of the last 
two centuries {2] . Among the various possibilities there is that of expressing the nth deriva- 
tive of a composite functions in terms of the Bell polynomials [6] , which are a generalization 
of the Hermite-Kampe de Feriet polynomials. A family of polynomials intermediate between 
that exploited so far and the Bell one is represented by the three variable Hermite-Kampe 
de Feriet polynomials, defined as Q 

The following identity generalizes Eq. fll.7p (see Refs. js] and 

D^e'^^' =i/(3n3ax2,3ax,a)e'^^'. (III.2) 
and, accordingly, it is evident that 

blf{x^) = Hf^^^^x" 4>.'i^4>A)^''"^ . (111.3) 

i.e., 

bif{x') = \n\Y Y ^^^^ ^^^V''"'" > /(x^). (in.4) 

I fc=0 m=0 ^ ' I 

It is evident that the method we are developing is a kind of modular scheme, which can be 
easily generalized and automatized. An example is provided by the following Leibniz rule 

bl\g{x'')h{x% = J2 (?) Hi;l,{Sx'^,3x%^)H^^\3x'fj,3xfj,fj), (III.5) 

k=0 ^ ^ 

where the addition formula 

Hi^\xi + X2,yi + y2,zi + z2) = (^) i^f2,(xi, yi, zi) iff (x2, ^2, ^22) (III.6) 

k=0 ^ ^ 

has been used. 
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IV. CONCLUDING REMARKS 



In this section we complete the previous discussion by using the properties of the higher- 
order Hermite polynomials. In addition, we will comment on other formulae reported in 



Ref. [l| about successive derivatives of functions. 

By using the identity Bar 



with 

[n/m] , 

we obtain the following generalization of Eq. ( 1III.3I) : 

[n/m] k2 N fc„_2 N 4 

Kl=0 k2=0 km-2=0 km-l=0 

»n— 1 m—1 

— In 

^ p=j q=0 / 

As stated in the incipit, this paper has been inspired by Eq. (1.1.1.1) of Ref. (l|. 
pass now to consider Eq. (1.1.1.2), which states that, for n > 1, 

4"/(y5)-E(-l)'ll^;^7^^:^'/(v/i). (IV,4) 



^ »n— 1 m—1 \ 



We 



-l)!(2Viy 

We can write 

V-4"/(v^) = e*^^/(y^) = /(v/^)=e^<^/o 



n! 

n=0 



and, taking into account t 
(see formula (25) of Ref. 



= e^^<^exp|-v^ + 0| (IV.5) 

le expression of the generating function for the Bessel polynomials 

3) 



n=0 k=0 ^ ' 



(IV.6) 



we obtain 



" . ^ 1 / t 



E ^^^/(v^) = E ^ ^"-1 ^ -^'/o. (IV.7) 
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from which Eq. (1IV.4P is easily obtained equating the coefficients of the same powers in t. 
As a final example of application of the method, we prove formula (1.1.1.3) of Ref. jl|, 

n / \ k—2n 

fc=0 ^ ^ )' 

Indeed, after setting 

X 



we can write 

4.n j-n 



oo 



n=0 n=0 



where in the second line we have used the identity [1 



e"^ ^-/(a;) = /(^^-^J , (|Ax|<l). (IV.IO) 

By remembering that [3] 

f;t"L„(x,y) = exp , (IV.ll) 



n=0 

where 



L„(..y)=.!g(-lf(;^^^ (IV.12) 



are the two- variable Laguerre polynomials, it's easy to show that (with L_i{x,y) = 0) 



n=0 



n \k J (n — k — 1)\ 

n=l k=0 ^ / ^ ^ 



from which, going back to the variable x, taking into account Eq. (11. 6|) . and comparing 
powers of t of the same order in OIV.Qp . Eq. OIV.Sp is obtained. 



Appendix A 

By setting x"^ = y in Eq. 01. 1|) we get the following operatorial identity 

[n/2] 



{V^D.y = J2 Si'^'\n, k) {^y-^' b--' (A.l) 



fc=0 
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where 



4^/2). ^ 7^- (A.2) 

^ ^ ' ^ 4''k\{n-2k)\ ^ ' 

Analogously, by setting a/x = z in Eq. 0IV.4p we get (n > 1) 

y-^^) =Y.^Si-'\n,k) z--~' (A.3) 



where 



The numbers S^\n^k) are a generalization of the Stirling numbers of second kind 
S2\n,k) [lot involved in the expansion 

n 

[xID,}j' = ^S'i\n,k)x''Dl. (A.5) 

A;=0 

Appendix B 

It is well known that symbols in mathematics have their own life and meaning. We can 
therefore generalize the operator defined in Eq. ( 111.81) as follows 

(2)f(")(a,6) = if(2)(2aa;7,67) (B.l) 

and note that 

n-2k l± \ 



If the function / can be expanded as in Eq. fll.2p . we can consider the problem of evaluating 
the following integral 



/oo 
dx(2)f(")(a/2,6)/(-c: 
-oo 

= (^J^ dxHi'\ax^,b^)e-'^^'^^ /o. (B.3) 

By remembering that the generating function of the two- variable Hermite-Kampe de Feriet 
polynomials is given by 

E^^^ni/,^)=e^'+^^ (B.4) 



ni 

n=0 
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we get 



°° j.n / poo 



/o 



(B.5) 



i.e., treating the integral on the r. h. s. as an ordinary Gaussian integral 



— X„(a, 6, c) = . exp 

n=0 » ' 



/o. 



(B.6) 



In this equation, by expanding the exponential in series and equating the t-like power terms, 
we find {n < 1) 

T2„(a,6,c) = 2"(2n-l)!!y^ (^1^ + 6^ /„_i/2, Xan+iK 6, c) = 0. (B.7) 

As an application of this result, from Eq. flll.QP one has 



dx (2)f(2")(«/2, &) ^o(a;) = 2"+^ (2n - 1)!! ^ 



T{n + 1/2)' 



(B.8) 
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